Tomographic characterization of three-qubit pure states with only two-qubit detectors 
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A tomographic process for tliree-qubit pure states using only pairwise detections is presented. 
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The understanding of multipartite entanglement is one 
of the objectives of the quantum information community. 
Bipartite entanglement is well understood for pure states, 
where the Schmidt decompositionQl plays the central 
role. Also much has been learnt about mixed states in re- 
cent years . The problem one has to face when the num- 
bers of parties grow is the consequent increasing in the 
complexity of the correlations involved. Recently, a sig- 
nificant progress in three-qubit correlations has been re- 
ported by N. Linden, S. Popcscu and W. K. WoottersQ, 
who have shown that ^"Almost every pure state of three 
qubits is completely determined by its two-particle reduced 
density matrices" . It is important to emphasize the con- 
text in which the work is developed which is understand- 
ing how information is stored in multipartite systems. 
Specifically, they show that there is no more information 
on the three-party state than what is already contained 
in the three reduced pair states. Their result is, at first 
sight, surprising. In the case of two qubits it is not dif- 
ficult to show that several global states give the same 
reduced states. For bipartite pure states, the Schmidt 
decomposition theorem0| asserts that, given a vector 
l^*) € V(8) W, one can choose orthonormal basis for 
V, and {\wj)} for W such that |*) = J2k \^k) ® \wk)- 
In the case of two qubits, V and W have dimension 2 and 
the vector state can be written asQ 

|^'(6i,(^)) ^cose\vi)(E)\wi)-he''^sme\v2)<E)\w2), (1) 

with 9 e [O, f ] and ip € [0, 27r]. Writing the density 
matrix in the basis {\vi) \wj)} makes it clear that the 
relative phase e"^ is locally inaccessible, i.e., for fixed 
9 all reduced density matrices are equal. So, there is 
more information on the full two-qubit state than on the 
parts represented by the reduced states. As a clarifying 
example, the two Bell statesH J$±) = {|00) ± |11)} /\/2 
originate the same local states [lOj. 

Suppose now that an experimental physicist wants to 
make tomographic measurements of three qubits with 
only two detectors. Reference 0] shows that all neces- 
sary information is available in the two-qubit reduced 
matrices, but does not suggest any practical method to 
do it. In this work we present a tomographic protocol for 



the complete characterization of generic three-qubit pure 
states, based only on pairwise detections. The protocol 
works whenever the result of Ref. Q holds, i.e.: for all 
pure states except for a restricted class that will be com- 
ment latter (see eq. It must be said that although 
coincidence measurements are allowed, only local opera- 
tions should be done, i.e.: one can think of three distinct 
laboratories, two of them cquiped with detectors at each 
time, and a classical electronic coincidence line between 
them. In this sense, we are only implementing LOCC 
(Local Operations with Classical Comunnication) . 

First of all, let us review the authors argument. Con- 
sider an arbitrary pure state = J2ijk ^ijh {"i-jk) of three 
qubits A, B, and C. A general state (pure or mixed) 
that has the same reduced states of can be obtained 
through a pure state j^"), describing three qubits plus 
an environment E (this process is called a purification, 
and its existence can be shown by the Schmidt decom- 
position). The fact that I^E") has the same reduced states 
than I J/) when restricted to each two-qubit subspaces puts 
restrictions in its form. It is shown that for a generic 
state 1^), these restrictions determine the form of j^*) as 
l^") = X^ijfc ^yfc Nj^)®!^); where the environment is fac- 
torized. Consequently the environment state is pure, and 
the three-qubit state is necessarily It is a very ele- 
gant argument, which also allowed Linden and Wootters 
to generalize this results for N qubits|^. It is, however, 
rather abstract, and gives no clue for the experimentalist 
to completely characterize his three-qubit pure state. 

The route we will take uses the generalization of the 
well known expression for one qubit: 
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where ctq is the 2x2 identity matrix, and ai are the Pauli 
matrices 
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where we leave blank all null entries, 
be obtained from the expression 

= Tr{po-^} . 



The coefficients can 



(4) 
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It is important to observe that eq. Q implies bo = I 
(normalization of p) and that the complete characteriza- 
tion of the state can be achieved with three mean value 
measurements b^, = (a^) (i.e.: by the three components 
of the so called Bloch vector). In fact, this is a tomo- 
graphic scheme for determining a spin-i state @. 
To generalize eq. (0} for three qubits, let us define 

S^^iy — (T^ ^ cTfj, ® a„, (5) 

and denote the state of three qubits by 

where we have adopted the convention of summation over 
repeated indexes throghout the paper (latin indexes from 
1 to 3; greek indexes from to 3). Once again, the co- 
efficients a^^j, can be obtained tomographically by the 
relation 

a^f^^ = Tr{pS^f^^} . (7) 

A first important consequence of eq. lUJ) is that aooo = 
1. As Pauli matrices are traceless, the reduced density 
operators are given by 

Pbc = Tr^ (pABc) = ^a-otiu^iiu, (8a) 
Pac = T^rsiPABc) ~ -^o-yQ^S^t,, (8b) 
Pab = Trc(pAsc) = 707^08-,^, (8c) 



where S^^ = cr,, ® a^. 

To directly determine p through eq. ((TJ one needs to 
evaluate sixty three mean values. Nine of them (3 x 3) 
are the three components of each Bloch vector (aioo, aojo, 
ooofc), and can be determined by individual detections. 
Twenty seven (3 x 9) are the pair correlations (0^0, a-Mk, 
aojk) and must be obtained through two-qubit coinci- 
dence measurements. The remaining twenty seven are 
three-qubit correlations, and are directly available only 
through three-qubit coincidence detections. For a general 
mixed state the number of mean values to be determined 
is exactly the same as the number of coefficients in the 
density operator. However, any previous knowledge on 
the state of the system is expected to reduce the number 
of parameters needed. In particular, for pure states we 
can use the idempotency relation, 

= P, (9) 

to obtain the coefficients a^fc from those available in the 
pair states. From expression ^ we get 64 equations, 
which can be organized in four sets: the first set consists 
of one solely equation 

J2^"''ioo+aojO+alok+a^jo+aMk+aojk+alk) = 7, (10a) 

ijk 

the second set is given by 

Sa.joo = a.ijoaojo + aiok.aook + a.ijkaojk, (10b) 
with similar equations under permutations of indexes. 



The third set is consituted of 

iO-ijO = OiOOaOjO + amkO-ijk + O-OikO-iOk — -^^ilt£jmu0.lmO(^tuO — -^^ilt£jrnu0.tukCllmk, (lOc) 

together with analogous equations under cyclic permutations, and where we use the Levi-Civitta symbol e^fc for the 
totally antisymmetric tensor. Finally, the fourth group 

(lOd) 



The tomographic process is thus constituted by the 
thirty six mean values measured in individual and two- 
qubit coincidences, and the sixty four equations (|10|l . 
that must be solved for aijk- The Linden, Popescu and 
Wootters' result guarantee the generic solution of the 
whole set of equations. 

Anyhow, the last set (jlOdp gives 27 linear equations 
on the 27 unknowns a^-fc. In case they are linearly inde- 
pendent, this specific set can give the complete solution. 
We numerically checked such independence in all of hun- 
dreds of random choices. However, as pointed out by the 
authors of Ref. Q, there are exceptions, for states like 

\GHZ{e,(p)) = cos 6*1000) +e*'^ sine* I 111) , (11) 



with 6 e (0, |-), in which the phase e"^ is pairwise inac- 
cessible, in the same sense as its analog in the two-qubit 
states QJ. Thus, all the exceptions are the states that, 
for some choice of local basis, can be written as (|ll|l . since 
for any other state, all the phases can be obtained with- 
out involving triorthogonal basis vectors. In such case, 
eqs. (|10d|l can not be linearly independent, but we con- 
jecture that a more geometrical argument can show the 
generic independence, and also point out the exceptions 

A parameter counting shows how rare the exceptions 
are. A vector in (g) ® is given by 8 complex num- 
bers (i.e.: 16 real numbers). Normalization and global 
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phase reduce it to 14 real numbers. Local unitary oper- 
ations are given by the action of SU (2) group in each 
aubit[lll|. SU (2) is parametrized by 3 real numbers (the 
3 Euler angles, or the 3 components of a vector v for the 
Lie algebraic paramctrization U (v) = exp {iv ■ a}) , so 
the orbit space has real dimension 5 (i.e.: 14 — 3 x 3)0 ■ 
As the fasc (p in (|11|) can be changed by local unitarities, 
the GHZ family has just one parameter, 6. So it is like a 
regular curve in a five dimensional manifold. 

The protocol here presented can be described in the 
following steps: first two-qubit measurements determine 
the 36 mean values characterizing all the two-qubit re- 
duced states {i.e.: the coefficients Oioo: ^ojo, aoofcj OyOj 
fliofc, and oojfc); then this experimental data is used as in- 
put on the 27 linear equations (|10d|l . and their solution 
gcnerically determines the 27 remaining coefficients a^fe, 
i.e.: generic three-qubit pure states are completely char- 
acterized by these 27 calculated coefficients plus the 36 
directly measured ones. In fact, we can check for the pu- 
rity of the measured state using the obtained coefficients 
to test the remaining 37 equations (|10alllOblllUc|l . If any 
of these "testing" equations is not satisfied (within exper- 
imental precision), one should conclude that the original 
state is mixed, and can not be determined with only 36 
mean values. 

Some other questions can be raised on this issue. Is 
there any other tomographic process, restricted to two- 
qubit detections, that can determine the state with fewer 
measurements, without introducing new exceptions? Is 



the optimal number of measurements, 14, achievable with 
this kind of restriction? Recently Diosi pointed out 
that a generic tripartite pure state can be determined by 
the knowledge of any two constituent pairs. Again this is 
a generic result in which interesting exceptions arise. For 
example, for three qubits, if an experimentalist decides 
to directly access pab and psc ■ and the prepared state is 
\^) = -^(1000) + |010) 4- e''^ |111)), it will be impossible 
to determine the phase ip. It would be available, however, 
at pac- In fact, it is an interesting problem to classify 
the multipartite pure states by the partial information 
necessary to completely determine them. Such a classi- 
fication could help understanding the curious geometric 
structure behind pure states. 

In this paper we provide for a feasible experimental 
prescription to completely characterize a generic three- 
qubit pure state using only two-qubit detectors. 
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